Highly anisotropic anomaly in the dispersion of the copper-oxygen bond-bending phonon in 
superconducting YBa 2 Cu 3 7 from inelastic neutron scattering 
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Motivated by predictions of a substantial contribution of the "buckling" vibration of the Cu0 2 layers to 
d-wave superconductivity in the cuprates, we have performed an inelastic neutron scattering study of this 
phonon in an array of untwinned crystals of YBa 2 Cu 3 7 . The data reveal a pronounced softening of the 
phonon at the in-plane wave vector q = (0, 0.3) upon cooling below ~ 105 K, but no corresponding anomaly 
at q = (0.3,0). Based on the observed in-plane anisotropy, we argue that the electron-phonon interaction 
responsible for this anomaly supports an electronic instability associated with a uniaxial charge-density 
modulation and does not mediate d-wave superconductivity. 



PACS numbers: 74.25.Kc, 74.72.Gh, 63.20.kd, 71.27.+a 

Research on the mechanism of high-temperature super- 
conductivity in the cuprates has recently made substantial 
progress based on a quantitative analysis of possible Cooper 
pairing interactions, which include coupling to spin fluctua- 
tions and phonons. The spin fluctuation mediated pairing in- 
teraction favors the experimentally observed d-wave pairing 
state, [lj but its strength relative to phonon-mediated pair- 
ing interactions remains a matter of intense debate. [fH 412ll 
Whereas density functional theory (DFT) indicates a small 
[00 or negligible 11101 lll || contribution of the electron- 
phonon interaction (EPI) to d-wave superconductivity, it has 
been argued that electron correlations can greatly enhance 
the coupling strength Q21 (3j such that it yields a substantial 
[1711 or even dominant [4] contribution to the pairing inter- 
action. In this context, much attention has been focused 
on phonons that modulate the Cu-0 bonds in the Cu0 2 
layers, which are generic to all cuprate superconductors. 
Most model calculations indicate that the EPI of high-energy 
stretching vibrations of these bonds is either detrimental 
[H or indifferent [JSHZl to d-wave pairing. However, 
several theories indicate that coupling to the lower-energy 
"buckling" vibration that modulates the Cu-O-Cu bond angle 
is attractive in the d-wave channel. ]5}j9l 

Two experimental strategies have been employed to test 
these theories and to quantify the strength of the different 
electron-boson coupling channels. The first approach has 
identified fingerprints of coupling to bosonic modes in the 
electronic spectral functions extracted from angle-resolved 
photoemission (ARPES), optical, and tunneling data, [TJO 
but these data are not specific enough to unequivocally 
discriminate between contributions from spin excitations 
and phonons. A second, complementary approach there- 
fore targets corresponding anomalies in the bosonic exci- 



tations. The prominent spectral-weight redistribution of 
spin excitations at the superconducting transition temper- 
ature determined by neutron scattering II 1 3 II indicates a 
positive contribution to the d-wave pairing susceptibility, 
in agreement with theoretical predictions. [1] The Cu-0 
bond-stretching vibration also exhibits anomalies indicative 
of a substantial EPI, H14H20II but most model calculations 
agree that coupling to this mode does not support d-wave 
pairing. These, however, have not predicted or reproduced 
the above-mentioned anomalies | |21 || (with the exception of 
Ref.lH). 

Much less information is available on the buckling mode, 
which has been predicted to favor d-wave superconduc- 
tivity. [JSHSI Raman scattering experiments on optimally 
doped cuprates have shown that the mode exhibits a 
superconductivity-induced softening of ~ 1.5 % at wave 
vector q = 0. 1221 12311 However, q-dependent data are 
required for a comprehensive assessment of the contribu- 
tion of the buckling phonon to the experimentally observed 
electron self-energy anomalies and to the d-wave pairing 
state. Such data can be collected by inelastic neutron scat- 
tering. Since the lattice dynamics is further complicated 
by structural disorder, prior neutron scattering work has 
focused on stoichiometric YBa 2 Cu 3 7 . H24II25II The results 
confirmed a modest softening of the buckling mode upon 
cooling below the superconducting transition temperature 
T c , which agrees with Raman scattering data 11221 12311 at 
q = and decreases rapidly with increasing q, in accord 
with model calculations. [26 , 27] However, a broadening of 
the neutron profiles with a maximum amplitude at nonzero 
q II25II was not explained by these calculations. Since these 
experiments were carried out on twinned specimens, the 
in-plane anisotropy of these features could not be resolved. 
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FIG. 1: (a-e) Constant-Q scans at Q = (0, 0.3, 10). The profile at energy 40-42 meV arises from the bucking mode. Panel (a) shows raw 
data. A sloping background (dotted line) was subtracted in the other profiles for clarity. The black and green lines are results of fits to a 
phenomenological expression comprising two Voigt functions and to a full resolution calculation, respectively, as described in the text, (f) 
Temperature dependence of the lineshape parameterized by the intensity of the two peaks in the phenomenological expression. 



In the light of recent reports of pronounced, strongly 
temperature and doping dependent in-plane anisotropies 
of the magnetic 1 1311 and transport 11281 |29ll properties of 
underdoped YBa 2 Cu 3 6+5 , we have used inelastic nuclear 
neutron scattering to determine the in-plane anisotropy of 
the dispersion of the buckling mode in an untwinned speci- 
men of YBa 2 Cu 3 7 . Our sample consists of 185 individually 
detwinned single crystals with total mass 2.63 g, co-aligned 
on an aluminum plate. Magnetization measurements on a 
representative set of crystals showed that T c = 90.0 ± 0.5 
K. X-ray powder diffraction measurements yielded the room 
temperature lattice parameters a = 3.8173, b = 3.8844, 
and c = 11.681 A. These data imply that the sample is 
nearly stoichiometric (5 ~ 1) and slightly overdoped. This 
is important because even a small density of oxygen va- 
cancies leads to the formation of superstructures that may 
affect the phonon dispersions. [30 ] The detwinning ratio of 
the array is 86 percent. The neutron measurements were 
performed on the PUMA spectrometer at the FRM-II reactor 
in Garching, Germany, and on the IT spectrometer at the 
ORPHEE reactor at LLB, Saclay France. The (2,2,0) reflec- 
tion of Cu and the (0,0,2) reflection of pyrolytic graphite 
were used to monochromate and analyze the neutron beam, 
respectively The final neutron energy was fixed at 8 mey 
and the data were taken in the Brillouin zones adjacent to 
the (1,0,10), (0,1,10), and (0,3,1) reciprocal lattice vectors. 
(Wave vectors are quoted in units of 2n:/(a, b, c).) 

Figure 1 shows representative neutron scattering pro- 



files in a scattering geometry that maximizes the structure 
factor of the buckling mode, at the in-plane wave vector 
q = (0, 0.3) where the broadening was observed in the early 
experiment II25II . In addition to the buckling mode at en- 
ergy E = 40 — 42 mey the data show a mode at 44.5 meV 
that is not strongly temperature dependent and will not 
be discussed further. The profiles of the buckling mode, 
on the other hand, exhibit a highly unusual behavior. At 
temperature T = 120 K, the profile shows a single peak 
with a barely discernible shoulder on the low-energy side 
(Fig. le). Upon cooling, however, the shoulder grows con- 
tinuously and draws intensity from the main peak, resulting 
in a two-wing profile at low temperatures (Figs. la-d). In 
order to parameterize the temperature evolution of the line- 
shape, we have fitted the data to a phenomenological form 
that comprises two Voigt functions with fixed energies 40.0 
and 41.8 mey and widths 1.6 and 1.2 meV (black lines in 
Figs. la-e). The fits yield a good description of the exper- 
imental data over the entire temperature range (although 
they do not correctly capture the physics; see below). The 
T-dependence of the resulting peak intensities (Fig. If) 
indicates a well defined onset of the lineshape anomaly at 
T = 105 ± 15 K. 

This unusual behavior is confined to in-plane wave vectors 
in the range q = (0, 0.2) — » (0, 0.4). The profiles for wave 
vectors outside this range along the same direction, and 
for the entire range along the perpendicular direction, q = 
(0, 0) — » (0.5, 0), are well described by a single, undistorted 
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FIG. 2: Constant-Q scans of the buckling mode at Q = (0.3, 0, 10) 
and (0, 0.3, 10). A sloping background (Fig. la) was subtracted. 



peak. A representative scan at q = (0.3, 0) is shown in Fig. 2. 
The strong in-plane anisotropy of the anomalous dispersion 
of the buckling mode revealed by these data could not be 
recognized in previous work on twinned crystals. 11241 12511 

Since the E — q range of the anomalous low-temperature 
dispersion of the buckling mode is comparable to the dimen- 
sions of the resolution ellipsoid of the neutron spectrometer, 
a careful resolution calculation is required to extract infor- 
mation from the observed lineshape. We have performed 
such a calculation based on the dispersion relation shown 
in Fig. 3, which is flat along the a-direction at all tempera- 
tures (Fig. 3a) and along b for T > 120 K (red line in Fig. 
3b), but develops a sharp dispersion anomaly along b at 
lower temperatures (Fig. 3b). Considering the simplicity of 
the model, the agreement between the calculated profiles 
(green lines in Figs. 1 b-e) and the data is satisfactory. In 
particular, the calculation provides an explanation of the 
two-wing profile discussed above. The low-energy wing 
results from the bottom of the dispersion anomaly at the 
nominal spectrometer setting. However, the instrumental 
resolution ellipsoid also encompasses a wide segment of the 
nearly flat dispersion surface away from the anomaly. This 
explains the high-energy wing of the profile. 

A highly anisotropic dispersion anomaly of the buckling 
mode that is strongly enhanced upon cooling below ~ 105 
K thus provides an adequate explanation of all aspects of 
the experimental data. In view of the presence of several 
other phonon branches in the energy range 35-45 meV 
previously detected by Raman and infrared spectroscopies 
II31II . one may wonder to what extent the T -dependent 
spectral-weight redistribution of the buckling mode results 
from mixing with other modes. We have calculated the neu- 
tron structure factors based on the eigenvectors provided by 
ab-initio DFT calculations II32II and found that most of the 
modes in this energy range (mostly in-plane vibrations of the 
apical oxygen ions, with some admixture of chain and plane 
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FIG. 3: Dispersion of the buckling mode along the a- and b-axes. 
The black line is the dispersion used for the resolution calculation. 
A projection of the four-dimensional resolution ellipsoid is shown 
for comparison. The data point at q = (0, 0.3) is the result of the 
resolution convolution described in the text, the remaining points 
were determined by fits to standard Voigt functions. The inset 
shows the eigenvector of the buckling mode at q = (0, 0.3). The 
elongations of the apical oxygen atoms and of the in-plane oxygen 
atoms along b were enlarged by a factor of four for clarity. 



oxygens) exhibit substantial spectral weight in the Brillouin 
zone around the reciprocal lattice vector (0, 3, 1). Scans in 
this zone for q = (0, 0.3) show no evidence of temperature 
dependent modifications of the phonon intensities or line- 
shapes outside the statistical error. While a spectral-weight 
transfer to other modes analogous to the one discussed for 
the Cu-0 stretching vibration II15II17II cannot be ruled out 
entirely, it therefore does not appear to play a major role 
for the buckling mode. This conclusion is also supported by 
the observation that the integrated intensity of the buckling 
mode profile shown in Fig. 1 is T -independent. 

We now discuss the microscopic origin of the dispersion 
anomaly. Since wave vectors around q ~ (0, 0.3) connect 
nested segments of the Fermi surface of YBa 2 Cu 3 7 , II33II 
an anomaly resulting from the EPI provides a natural expla- 
nation of our data. The strong temperature dependence of 
the anomaly (Fig. 1) can then be attributed to the opening 
of a gap (or pseudogap) on the Fermi surface. Although the 
data of Fig. If are consistent with a sharp onset or enhance- 
ment at T c , the momentum-space structure of the anomaly 
speaks against a direct relation to the superconducting gap. 
11231 |34l 13511 Raman scattering experiments have shown 
that the pair-breaking energy 2 A ^ 60 mey with a ~ 10 
% anisotropy between a- and b-directions. 11231 13411 Since 
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2A significantly exceeds the phonon energy along both a- 
and b-axes, the onset of superconductivity is expected to 
affect the phonon lineshapes similarly in both directions. 
Specifically, a softening of the phonon is expected as the gap 
sweeps through the phonon energy upon cooling. 11361 13711 
This is indeed the case for the superconductivity-induced 
softening centered at q = (Fig. 3), []22l [23] but not for 
the anomaly at q ^ 0, which is present only along the b-axis. 
The strong anisotropy in phonon softening is indicative of a 
uniaxial EPI, which cannot be responsible for the supercon- 
ducting gap whose magnitude is similar along a and b, as 
further explained below. 

Based on these considerations, a (real or incipient) in- 
stability other than superconductivity must be responsible 
for the T-dependent anomaly of the buckling mode. The 
strong in-plane anisotropy of the anomaly indicates that it 
has a one-dimensional structure, which may arise from a 
uniaxial charge density wave (CDW) instability. There is 
independent evidence for a CDW originating in the CuO 
chains in the YBa 2 Cu 3 7 crystal structure, tl3"8ti40l and the 
wave vector of the phonon anomaly is similar to the one 
of the CDW recently reported by x-ray diffraction. [40] 
The temperature onset of the phonon anomaly (Fig. If) 
is also near the onset of a pronounced in-plane anisotropy 
of the Nernst effect in YBa 2 Cu 3 7 , H29II which has been 
interpreted as evidence of a stripe or Pomeranchuk instabil- 
ity originating in the Cu0 2 planes. We note, however, that 
the wave vector of the instability indicated by the phonon 
anomaly is perpendicular to the one characterizing the uni- 
axial spin density wave state recently observed in strongly 
underdoped YBa 2 Cu 3 6+5 . 1141 II This may imply that spin 
(charge) driven instabilities with propagation vector along 
a (b) dominate the interplay with d-wave superconductiv- 
ity in the low (high) doping regimes of the YBa 2 Cu 3 5+5 
phase diagram. Such a multi-phase competition may be 
influenced by the hybridization between the chain-derived 
energy bands with those arising from the Cu0 2 layers. 1 1 3 3 II 

We finally address the magnitude of the dispersion 
anomaly. If we assume that the phonon softens due to 
an opening of a gap on the Fermi surface, (see Ref. [36] 
for the case of the superconducting gap, although our ar- 
gument would work for a gap of any origin) the phonon 
renormalization magnitude should be proportional to the 
q-dependent EPI strength. The observed ~ 6 % softening of 
the buckling phonon at q ~ (0, 0.3) upon cooling indicates 
that the EPI of the buckling mode may be strong enough to 
contribute significantly to the electronic self-energy along 
the antinodal direction of the d-wave gap function, as hy- 
pothesized based on ARPES measurements. []2l [3 Our data 
provide an improved basis for a quantitative assessment 
of this contribution relative to the one generated by spin 
fluctuations 11 3 II . Despite its strength, we have shown 
that the highly anisotropic EPI responsible for the anomaly 
at q ~ (0, 0.3) does not support d-wave superconductivity. 
We conclude that the buckling mode does not contribute 
to d-wave pairing beyond q = 0, where the phonon soften- 



ing below T c translates into a small dimensionless coupling 
constant X d ~ 0.02, in agreement with DFT calculations. 
| | 1 1 1| It is possible, however, that it contributes to the ob- 
served sub-dominant s-wave admixture to the pairing state. 

11251 1551 1551 

In any case, the sharp phonon anomaly indicates that 
overdoped YBa 2 Cu 3 7 is close to an instability associated 
with a uniaxial CDW. Further work is required to assess the 
contribution of the electron-phonon coupling driving this 
instability to the anomalous band dispersions revealed by 
ARPES. 
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